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Abstract 

We have considered the conductivity properties of a two dimensional elec- 
tron gas (2DEG) in two different kinds of inhomogeneous magnetic fields, i.e. 
a disordered distribution of magnetic flux vortices, and a periodic array of 
magnetic flux vortices. The work falls in two parts. In the first part we show 
how the phase shifts for an electron scattering on an isolated vortex, can be 
calculated analytically, and related to the transport properties through the 
differential cross section. In the second part we present numerical results for 
the Hall conductivity of the 2DEG in a periodic array of flux vortices found by 
exact diagonalization. We find characteristic spikes in the Hall conductance, 
when it is plotted against the filling fraction. It is argued that the spikes 
can be interpreted in terms of "topological charge" piling up across local and 

global gaps in the energy spectrum. 

PACS 73.50.-h, 73.40.Hm, 73.20.Dx, 73.50.Jt 
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I. INTRODUCTION 



Over the last decade the two dimensional electron gas (2DEG) have been exposed to a 
wide range of physical experiments, in which the electrons have been perturbed by differ- 
ent configurations of electrostatic potentials, with or without a homogeneous perpendicular 
magnetic field. These experiments have shown new kinds of oscillations in the magnetocon- 
ductivity, with a periodicity not given by the geometry of the Fermi surface, as is the case 
with the Shubnikov-de Haas oscillations, but given by the interaction of the two length scales 
given respectively by the magnetic length, and by the spatial structure of the potential, e.g. 
the Weiss oscillations jT|. More recently, there have been increasing interest in systems where 
the 2DEG is exposed to an inhomogeneous perpendicular magnetic field. In such systems 
the inhomogeneities in the magnetic field acts as perturbations of the 2DEG, relative to 
the homogeneous magnetic field, where the band structure consists of the completely flat 
Landau bands. The inhomogeneous magnetic field appears in the Hamiltonian in the form 
of a non trivial vector potential. In the case of a periodic variation in the magnetic field, it 
is possible to construct a periodic vector potential, if and only if the flux through the unit 
cell of the field is equal to a rational number, when measured in units of the flux quantum 
0o = h/e. Under these special circumstances the Hamiltonian is periodic, and Bloch states 
can be used as a basis for the calculation of responce properties of the electron gas. 

In this paper we have considered a special class of spatially varying magnetic fields 
which consists of flux vortices, that are either distributed at random or placed in a regular 
lattice structure. A system consisting of a 2DEG penetrated by a random distribution of 
magnetic flux vortices, have been experimentally realized by Geim et al. H||. They made 
a sandwich construction of a GaAs/GaAlAs sample with a 2DEG at the interface, and 
a type II superconducting lead film, electrically disconnected from the 2DEG. When the 
system was placed in an external magnetic field, and cooled below the transition temperature 
of the film, the magnetic field penetrated the film, and thereby the 2DEG, in the form of 
Abrikosov vortices. When the external magnetic field is weak, below 100G, the vortices will 
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be well separated, and the 2DEG therefore sees a very inhomogeneous magnetic field. In the 
experiments conducted by Geim et al. the flux pinning in the film was strong, resulting in a 
disordered distribution of flux vortices. This is the physical situation which we investigate in 
Sec. [Tl] below. In very clean films of type II superconducting material, the flux vortices will 
order in a periodic array, i.e. an Abrikosov lattice, and thereby create a periodic magnetic 
field at the 2DEG. This is the situation which we analyse in Sec. |TTT| . 

Several authors have investigated the transport properties of 2DEG's in different kinds 
of inhomogeneous magnetic fields. Peeters and Vasilopoulos Q have made a theoretical 
study of the magnetoconductivity in a 2DEG in the presence of a magnetic field, which was 
modulated weakly and periodically along one direction. They found large oscillations in the 
longitudinal resistivity as a function of the applied magnetic field strength. These oscillations 
are due to the interference between the two length scales given respectively by the period 
of the lateral variation of the magnetic field, and by the magnetic length corresponding to 
the average background field. The oscillations are reminicent of the Weiss oscillations, but 
have a higher amplitude and a shifted phase, relative to the magnetoresistance oscillations 
induced by the periodic electrostatic potential. 

The problem of how the transport properties of the 2DEG is modified by the presence 
of a random distribution of flux vortices, have been treated earlier by A. V. Khaetskii ||, 
and also by Brey and Fertig |J. The approach used by these authors are basically similar to 
the one we have presented in Sec. |I|, i.e. based on the Boltzmann transport equation. The 
main difference being that while Khaetskii have treated the scattering in certain limiting 
cases, including the semiclassical, and Brey and Fertig have calculated the scattering cross 
section numerically, we have found an analytic expression for the scattering cross section for 
electrons scattering on an idealized vortex. This has enabled us to study the scattering in 
more detail, and to observe scattering resonances. 

In Sec. [TTT] we will address the "paradox" of how the Hall effect can disappear in the 
following situation: We imagine a 2DEG in a regular 2D-lattice of flux vortices, with the 
magnetic field from a single vortex exponentially damped with an exponential length £, 
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in units of the lattice spacing. We take the total flux from a single vortex to be 4>o/2, 
as is the case when the vortices come from a superconductor. When £ ^> 1, the field 

2 

is homogeneous and the Hall conductivity is ou = Pjr, where p is the number of filled 
bands. In the other limit i.e. when £ — > 0, the time reversal symmetry is restored, and the 
Hall conductivity vanish. The fact that the system has time reversal symmetry when the 
vortices are infinitely thin, can be seen by subtracting a Dirac string carying one quantum 
of magnetic flux <f>Q, from each flux vortex. The introduction of the Dirac strings can not 
change any physical quantities, and the procedure therefore establishes that the system, 
with infinitely thin vortices with flux +0o/2, is equivalent with the system with reversed 
flux — 0o/2, through each vortex. The paradoxical situation arises because it is known from 
general arguments JF|, that the contribution to the total Hall conductivity from a single filled 
nondegenerate band is a topological invariant, and therefore cannot change gradually. The 
situation is even more clear cut if we imagine a periodic array of exponential flux vortices 
carying one flux quantum O each. Then the limit £ = oo corresponds to a homogeneous 
magnetic field, while the opposite limit £ = coresponds to free particles. Incidentially this 
scheme can be used to establish an interpolating path between the multifractal structure 
known as Hofstadters butterfly || which is a plot of the allowed energy levels for electrons 
on a lattice in a homogeneous magnetic field, and the corresponding plot for lattice-electrons 
in no field, which is completely smooth ||. 

The plan of the paper is as follows. In Sec. |H] we concentrate on the theory of electrons 
scattering on a single vortex, and the physical consequences for the resistivities. First we 
review the classical scattering theory in Sec. [II A| , before we discuss thre quantum theory of 
scattering in Sec. [II B| . The longitudinal and transverse resistivities are discussed in Sec. |I C| , 
and resonance scattering is demonstrated in Sec. |11 D. The case of a 2DEG in a periodic 
array of flux vortices, is the subject of Sec. |T|. In the first part, Sec. [Ill A| -P flTJ| , of this 
section the general theory of electron motion in a periodic magnetic field is reviewed, and 
in the second part, Sec. [ill h^ITiT] , we present the results of the numerical calculations. 
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II. SINGLE VORTEX SCATTERING 



In this section we will consider the consequences of the introduction of magnetic flux 
vortices into a 2 dimensional electron gas, in the approximation where each vortex is treated 
as an individual scattering center. The vortices are assumed to be distributed at random, 
homogeneously over the sample. The average separation between the vortices is assumed so 
large, that we can neglect interference from multiple scattering events. In the experimental 
situation the mean free path If from impurity and phonon scattering, may be very long 
compared to the average separation between the vortices, due to the very clean samples and 
liquid Helium temperature. This means that multiple scattering and interference may have 
important consequences. Nevertheless we will stick to the simplifying picture of vortices 
as individual scatterers in this section. As we shall see, the gross features observed in 
experiments on this system, can be accounted for within this approximation. 



A. The Classical Cross Section 

We will start by calculating the differential cross section for an electron scattering on a 
flux vortex within the framework of classical mechanics. This will provide a reference frame, 
and allow us to speak unambiguously about the classical limit. In the calculations we shall 
use an ideal vortex, which has a circular cross section with constant magnetic field inside, 
and zero magnetic field outside 



B(r) = < 



B o = for r < R v 

-It- (x) 



"V 

for r > R„. 



Here R v is the radius, and is the total flux carried by the vortex. The classical orbit is 
found as the solution to Newton's equation of motion with the force given by the Lorentz 
expression F = —ev X B. It consists, as is well known, of straight line segments outside the 
vortex, and an arc of a circle inside, with radius of curvature given by the cyclotron radius 
l c = v/cu c , with v being the particle velocity, and u c = — the cyclotron frequency. Inside 



the vortex the orbit is an arc of a circle, and as it is impossible to draw a circle that only cut 
the circumference of the vortex once, it is a simple geometrical consequence that a particle 
obeying the laws of classical mechanics, and which initially is outside the vortex, can never 
become trapped inside the vortex. It is clear that the classical scattering is controlled by 
the single parameter 7 = l c /R v , which is the ratio between the radius of the cyclotron orbit, 
and the radius of the flux vortex. Our first objective is therefore, for a given 7, to find 
the relation between the impact parameter b, and the scattering angle 9. Fig. [I] shows the 
geometry of the scattering, and the definition of the impact parameter b, and angles 0, ijj, 9. 
Let us define the reduced impact parameter j3 = b/R v , which is bounded to the interval 
— 1 < /3 < 1. By inspecting Fig. [I] it is observed that the following relations hold 

P = sm(p (2) 

Q 

7 sin - = sin(V> - 0)sign(7 + (3), (4) 

where the sign of 7 is dictated by the direction of the magnetic field inside the vortex. We 
take 7 to be positive. After a small amount of arithmetic and ip are eliminated, and we 
have 



*4 = sign(7 W 7 2 + w+i - (5) 

This relation gives the scattering angle as a function of the impact parameter. In classical 
scattering the scattering angle is always uniquely determined, once the impact parameter is 
given, in contrast to the inverse, i.e. the impact parameter as a function of the deflection 
angle. 

In an experiment one would measure the number of particles per time — j L -^-d9 scattered 

dt 

to an interval d9 about the angle 9. This will of course depend on the incoming flux of 
particles j, defined as the number of particles per time that cross a unit length perpendicular 
to the flow. Therefore we write 

dN{9) .da 
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The differential cross section 4| gives the total weight of impact parameters, which give 
scattering into the direction 0. If, for a given angle 9, we label the different values of the 
impact parameter b\, b%, . . . b p , which result in scattering into 9, then we have 



da 
d9 



E 

i=l 



d9 



(7) 



Equation R has at most two solutions, which are easily found to be 



/3±(6) = -7 sin 2 0/2 ± cos 0/2 /l. - 7 2 sin 2 9/2 



(8) 



The solutions have to obey the auxiliary conditions \(3\ < 1, and sign(7 + (3) =sign(0). 
Furthermore we have 



d(3 d 



-7 cos 0/2 sin 9/2 =f sin 0/2- 



1 + 7 2 cos 



(9) 



dd 2y/l -7 2 sin 2 0/2 
from which the differential cross-section can be calculated from Eqn. [7], still having the 
auxiliary conditions in mind. Examples of cross-sections and trajectories are shown in 
Fig. 0. The integrated cross section 



0"tot 



f"K 


da 


\ d9 


I —n 


d9 



(10) 



is equal to the total weight of impact parameters, which hit the vortex. It is equal to 
the diameter of the vortex a tot = 2R V , as is always the case in classical scattering. Let 
us imagine an electron at the Fermi surface scattering off the vortex. Then we have l c = 
vp/io c = mvp/eB = hkp/eB. If we furthermore take the flux of the vortex to be a fraction 
/ of the flux quantum O , the flux density becomes B = (/0 o )/7ri2 2 = 2fh/eB 2 v . The 
dimenensionless cyclotron radius is then given by 7 = l c /R v = kpR v /2f, and this parameter 
we call n/2f. In the quantum regime k = 2ttR v /Xf, and / are the natural parameters 
to characterize the scattering. This identification of parameters allow us to compare the 
predictions of classical and quantum theory below. 
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B. Quantum Scattering 



In this section we shall consider the electron scattering off a magnetic flux vortex within 
the framework of quantum scattering theory. We will calculate the differential cross section, 
from which the longitudinal and transverse conductivities can be found from the theory of 
Sec. |11 Q The quantum nature of the electron radically alters the picture of the scattering 



process, when the wavelength of the electron is comparable to, or longer than the diameter 
of the vortex. In the limit of very small electron wavelength, the scattering can essentially 
by described by the laws of geometrical optics, and thereby classical mechanics. 

We will again take an idealized cylindrical vortex, with constant magnetic field inside, 
and zero field outside, Eqn. [p. This vortex is completely symmetric under any rotation 
about the center axis. This symmetry can also be made a symmetry of the Hamiltonian, by 
chosing a proper gauge when writing down the vector potential. In cylindrical coordinates 
A = e r A r + egAg, we have 

B{r) = d r A 9 - -d e A r + (11) 
When B is invariant under rotation, this equation has the simple solution 

A r = 0, A e {r) = - f dr'r'B(r'), (12) 



r Jo 



which in our case give the vector potential A = egAg with 

{(fir 
for r < R v 
2**2 (13) 
- — for r > Rv. 

The Hamiltonian is given by the expression 

Here we have taken the charge of the electron to be — e. The Hamiltonian is rotationally 
invariant, and therefore commutes with the angular momentum about the symmetry axis, 
L z . Consequently L z and H have common eigenstates. The canonical momentum of a 
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charge-g particle in a magnetic field is given by the expression p = mv + qA = |V, and 
the operator for the angular momentum about the z-axis is L z = [r x p] z = ^-d$. The 
eigenstates of L z are e de , and the requirement that the wave function have no cut when the 
vector potential is non-singular, reduces the possible values of I to the set of positive and 
negative integers. We can now separate the variables of the common eigenstates of L z and 
H, and write 

Mr,0) = Rki(r)e ae , (15) 

where k is an energy label E = Let us introduce the flux quantum <f) = h/e and 

the dimensionless fraction / = (f)/<po- The differential equations for the radial part of the 
wave function, takes a particularly simple form if we write it down in dimensionless variables 
£ = r/R v and k = kR v = 2irR v /X. The energy variable k measures the size of the vortex 
compared to the electron wavelength. In terms of n and / the classical limit will be k, f ^> 1. 
With these definitions, the equation for the radial part of the wave function for £ < 1 is 

R"+ 1 -R'+(^ 2 -( l - + 10^ = 0- (16) 

And for £ > 1 we have 



R" + i# + U 2 - )r = q. (17) 



Inside the vortex an analytical solution to the radial equation can be found by the following 



procedure essentially due to L. Page pO| , pT[ . First we make the substitutions p = v2/£ and 
id 

2./' 



2 

w = 2? ( we assume / > 0), which results in the equation 



«• + lg + („-(! + |) 2 ) R = , (18) 

for < p < y/2f '. Next we write the radial function as 

R^p) = p m e-o 2 l%{p), (19) 
with m= |Z|, and insert this into Eqn. [18]. Hereby we get an equation for Vi 



/ 2m 4- 1 \ 

V/'+i J - p\V{ + {w-l-m- l)Vi = 0. (20) 

This equation can be further simplified by making the substitution x = p 2 /2. Finally we 
have the equation 

xV{' + (m + 1 - x)V{ - -{m + 1 + 1- w)Vi = 0. (21) 

This differential equation belongs to a class of equations known as Rummer's equation. 
Rummer's equation is a member of an even bigger class of equations of the form J2p=oi a p + 
b p x)^ = 0, which can all be solved by Laplace's method |JI2] . Kummer's equation is solved 



by the confluent hypergeometric functions M and U (in the notation of Abramowitz and 
Stegun |13| ). The complete solution to Eqn. |2l] can be written 



Vi(x) = c 1 M^-{l + m + 1 - w),m+ l,x) + c 2 U(^(l + m + l-w),m+ l,x). (22) 

It turn out that in order that Ri((,) be regular as £ — > 0, we must take c% = 0, and we can 
therefore write down the solution to the Schrodinger equation inside the vortex 

mw) = c^e-y? M( l -{i +\i\ + i- ^), i/i + ijey w . (23) 

Here G\ is a normalization constant which we will not need to evaluate. 

Outside the vortex the radial equation is just the differential equation for ordinary Bessel 
functions of the first kind. We therefore immediately have for £ > 1 

9) = (A t J l+f (K$ + AWO) j W - ( 24 ) 

The two constants Ai, Bi are found from the requirement, that the wavefunction has to be 
continuously differentialble at the boundary of the vortex. There is no need to normalize 
the wave functions <ft K i, as the normalization constant will drop out of the final expression. 

In quantum scattering theory one seeks a particular eigenstate of the Hamiltonian which 
belongs to the continnuous part of the spectrum, and which far away from the scattering 
center has a direction in which it represents an incoming flow of particle current. Let us 
consider an eigenstate coresponding to the energy E k - 



2m 
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i 



(25) 



We want ipk to represent the scattering of particles which are incident along the x-axis, as 
indicated in Fig. |^. The asymptotic boundary condition on ip^ is therefore that it far to the 
left of the origin represents an uniform current of incoming particles. 

The vector potential gives a contribution to the current j, so that the plane wave in 
the direction of the x-axis is altered from the field free form e lkrcose _ The particle current 
density is given by 



h 

2mi 



n 



n 



t 



(26) 



and it is straightforward to check that the correct form for a state with uniform current in the 
direction of the x-axis is e tkrcose - t f 9 . if the flux through the vortex is not an integer number 
of flux quanta, i.e. if the fraction / is not integer, then the factor e -4 -' is not single valued 
as it stands, and we have to introduce a cut to make it so. We only enforce the asymptotic 
boundary condition along the negative x-axis, so the cut can be placed anywhere outside 
this region. Let us for a moment introduce the principal angle [9], defined by [9] = 9 when 
c < 9 < c + 2ir, and otherwise given by periodicity. Then the single valued factor e 1 ^ has 
a cut along the half line 9 = c. We want to express the plane wave, as a sum over partial 
waves, and therefore we consider the inner product with e , I integer 



d9 
— < 

2tt 27T 



,ikr cos 6 -if \G\-iW 



c + 2 ^ d9 , 



ikr cos e-i{f+t)6 



2vr 

c +¥ d9 , 



c+. 



2tt 



(27) 



where the last equality sign holds if c 



3tt 



. This means that the cut is placed along the 



positive y-axis, as shown in Fig. |3|. The function J u {z) is known as Anger's function |14|, and 
coincide with Bessel's J u (z) when v is an integer. For v not an integer the Anger function 
has the nice property that is goes asymptotically as the Bessel function for large arguments, 



i.e. 



J„(z) = J u {z) + 



sin ixv 



7TZ 



l-- + 0(\z\ 



(28) 
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for z i — > oo. Let us remark that if the cut is placed somewhere else the integral will not give 
Anger's function, but asymptotically it will still go as some combination of Bessel functions. 
Let us now subtract the plane wave from -0& 

Mr, 0) - e ikrcosd -^ e = -i l +fj l+f {Kt) + + biBft+fiKt). (29) 

In the asymptotic region the involved functions can be expanded to give 

Vfc(r,0)-e* Mfl ^ = — L= {IkA-i^f-ib^e^+m-^ + 

M, - i l+f + ibfr] e -«e+W)l + l f } . (30) 



This combination of terms are the sum of an incoming and an outgoing circular wave. The 
coefficient multiplying the incoming wave must vanish, as all the ingoing current should be 
represented by the plane wave. We therefore have the condition 

hAi - i l+f + ib x B x = 0. (31) 

The radial differential equations are real and linear, and therefore A\ and Bi are real numbers 
by construction. The phase shifts Si are defined by the relation 

Si = arctan — . (32) 
Ai 

We note that the phase shifts are independent of the arbitrary normalization of the wave 
functions. If we write Ai = Ci cos Si and Bi = Ci sin Si, then Eqn. [31] can be solved to give 
biCi = i l+f e- iSl . The out going circular wave which represents the scattered current, is given 
by the expression 

p ikr ikr 

Putting things together, we get the following expression 

Fi = -J— e^~ i5l smSi. (34) 

V 7TK 

It is seen that Ti is a function of I only through the phase shifts Si, we can therefore write 
Ti = J 7 [Si\. From the above expression we can in principle calculate ^F(&), and thereby the 
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differential cross section 4| = \F(9)\ 2 . But this is only practically feasible if the sum over I 
converges so fast that we can approximate it with a finite sum. The way things are stated 
above, T\ does not go to zero for / — > — oo, but rather goes to a constant value. The reason 
for this is that we have not singled out the Aharonov-Bohm contribution to the scattering 
amplitude T[Q\ which is of a singular nature. We will now pass to the Aharonov-Bohm 
limit, that is the limit where R v — > 0, while the flux is kept constant. We will then be able 
to express the effect of the finite radius of the vortex, as the difference in the scattering 
amplitude from the Aharonov-Bohm result 



When treating potential scattering for potentials with finite range, the scattering wave 
function is directly written as the sum of the incoming plane wave and the outgoing circular 
wave. The same thing can not be done here, as can be seen by the fact that such a sum must 
have a cut, while the true wave function can not have any cuts. This is a consequence of the 
long range of the vector potential which far away from the scattering center falls off as 1/r. 
The term F(8)e lkT ' j \fr must therefore be interpreted as the largest term in an asymptotic 
expansion of the part of the wave function carrying the outward particle current ||15|| . This 
is not different from the procedure used to calculate cross sections for scattering on long 
range scalar potentials, i.e. Coulomb. 



In the limit of vanishing radius of the flux vortex R v = 0, we are left with only one 
dimensionfull variable k, and therefore it is impossible to express the scattering amplitude 



F(0) = 5F(6)+F AB (9). 



(35) 



1. The Aharonov-Bohm Limit 



in dimensionless form. When R, 



V 



the eigenstates of the Hamiltonian are everywhere 



given by the expression 




(36) 
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The physical demand that we have to impose on the solutions 4>m is that of boundedness at 
the origin. Let us put v — I + f, the order of the Bessel functions. The properties of J u (z) 
and Y u (z) for z — > 0, we infer from the relations 

u>) = i?rx ,J;i> +1) (37) 

= 4M^lM , (38) 
sin vix 

valid for all values of v and z. For v > J v is regular, and Y"„ is irregular, and we therefore 
have to take A\ = 1, Bi = for / + / > 0. For v < the condition for regularity is seen to be 
B1/A1 = —tamru. In order to select one of the two solutions of tan5f 4jB = — tan7r(Z + /) = 
— tan 7r/, we impose the additional condition that the radial wave function shall be positive 
around r = 0. Which translates into (—1)' sin 5f B sin ixf < 0. All in all this means that the 
phase shifts in the Aharonov-Bohm limit are 

, D f for l + f > 
S? B = (39) 
I — tt(Z + /) for I + f < 

modulo 27T. We note that it is only the fractional part of / which play a role in this limit, 
any integer part may be absorbed in a redefinition of I. We also note that the phase shifts 
are independent of the electron wavelength. If we take < / < 1, the calculation of the 
scattering amplitude look like this 

F AB (6) = i/4^ +i7r/ sin7r/f>- ae = -^e-^f-^^f. (40) 
V ; V nk £i sin § V ; 

The differential cross section for scattering on an infinitely thin vortex carrying a magnetic 
flux f<p is therefore 



da 

de 



1 sin 2 ./ (4i) 



4B 2nk sin 2 9/2 

A result first derived by Y. Aharonov and D. Bohm in 1959, by a slightly different ap- 
proach [jnj . The Aharonov-Bohm cross section corresponding to / = 1/2 is plotted in Fig. |4] 
for reference. The AB cross section is completely symmetric under reflection 9 1— > —9, for 
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all values of /, and can not give rise to a net transverse force on the electrons. On the other 
hand in analogy with ordinary impurities, it can give rise to a finite lifetime of the electron 
states, and thereby give a contribution to the longitudinal resistance of the system. We note 
that the cross section is periodic in the flux, with period equal to the flux quantum, and 
that it completely vanish when the flux is equal to an integer number of flux quanta. This 
fact, that the scattering amplitude F AB {6) disappers when the stringlike vortex contains an 
integer number of flux quanta, has useful consequences. 

The AB cross section is non integrable because of the singularity at 9 = 0, and we can 
therefore not calculate a tot . This is due to the long range nature of the interaction, i.e. the 
vector potential, which fall of only as 1/r. This situation is not different from scattering on 
an ordinary scalar potential with long range, such as the Coulomb potential. 



2. Phase Shifts for Scattering on Vortex with Finite Radius 



In this section we will briefly describe how the phase shifts for scattering on a cylindrical 
vortex with finite radius are calculated. The equation from which the phase shifts are 
derived, is simply the equation which results from the demand that the logarithmic derivative 
of the radial wave function must be continuous at the boundary of the vortex 

1 dRf 



J_ dRf 

Rf^f 



(42) 



Inside the vortex we have 



F _ 1 dRf 



\A-f + V 



aiM( ai + l, k + lj) 
k M(a h kJ) : 



where we have defined the parameters ai — ^(l + \l\ + 1 — fj), and hi 
vortex the logarithmic derivative reads 



(43) 

1. Outside the 



1 dR> 



R> d£ 



ji + yi tan S t 



Ji+f( K ) + >H-/(K)tan<V 
where we have introduced the abbreviations ji = kJ/ + j_i(k) — (I + /)J/ + /(/t) and yi 
KYi + f-i(K) — (I + f)Yi + f(K). It is now simple to solve for Si 



(44) 
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, r Jl ~ ElJl +f (K) 

tan ft = — . r • 45 

EiYi +f (K) -yi 

The tan ft 's are the bricks from which cross sections and transport coefficients can be build 
up. The presented curves and cross sections have all been calculated with phase shifts found 
by this expression with the help of Mathematica, which have implementations of all the 
involved special functions. 



3. Asymmetric Scattering on Vortex with Finite Radius 

For a vortex with finite radius, we want to express the scattering amplitude T{9) exclu- 
sively in terms of the dimensionless parameters k and /. The procedure for calculating the 
scattering amplitude, for scattering on a vortex with finite radius, is as follows. First the 
SFis are found from 

5Ti^T[5 Kl ]-T[5t B l (46) 

where 5 K i is the phase shift calculated numerically by the formulas given in the above section, 
and is given by Eqn. |34|. For small values of k the ftF^'s vanishes rapidly when \l\ 
increases, making it possible to approximate the sum J2i SJ^ie 119 well by a finite number of 
terms. We have found that for k < 10, of order 20 terms are needed at most. The scattering 
amplitude is then simply given by 

/max 

F{0)= £ W 16 + F AB (9), (47) 

where the dimensionless AB cross section is given by JF AB (6 I ) = J2i •7 r [ft AB ]e^ 6 '. From which 
we get the dimensionless cross section 

^ = W)| 2 - (48) 

Plots of differential cross sections calculated by this procedure is shown in Fig. |5|. We note 
that the degree of asymmetry is determined by the size of the parameter k = kR v . The 
classical limit is approached when k ^> 1. The quantum cross section, unlike the classical 
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one, gives a finite probability of scattering to both sides of the vortex for all parameter 
values, but not necessarily in all directions. 



C. Conductivity of 2DEG in Vortex Field 

In this section we will make an estimate of the contribution to the resistance of a 2DEG, 
from a random distribution of flux vortices. The experimental situation we have in mind, 
is that of Geim et. al. who placed a thin film of lead on top of a GaAs/GaAlAs 

heterostructure. When the temperature is lowered below the critical temperature, and the 
magnetic field is below H c2 , the magnetic field penetrates the superconductor in the form 
of Abrikosov vortices. At the 2DEG the field will be confined to regions of radius (or 
exponential length) A s , each threaded by a magnetic flux 4>o/2. The vortices are assumed to 
be distributed randomly, with an average separation d given by the strength of the external 
.B-field according to the relationship d 2 B ext = 4>o/2. 

The Boltzmann equation, linearized in the external electric field, reads 

,df° f d 2 q f£ ... , f p -f 



- ev ■ E ~fc = J (2^)2 {fp+« W P+«-*P ~ fp w P^p+i} ~ T . mp ■ ( 49 ) 
Here the transition probabilities w p ^ p+q are potentially asymmetrical quantities, due to the 



time reversal breaking magnetic field in the vortices. As argued by B. I. Sturman the 
correct form of the collision integral, even in the absence of detailed balance, is the one given 
in Eqn. f|9[ The electron-vortex scattering will be elastic. In order to solve the Boltzmann 
equation we Fourier transform, and write 

oo 

/(M)= E ^ in6 fn(k) (50) 

n=— oo 
oo 

w(k,9)= E e~ me w n (k), (51) 



where 9 is the angle between k and E. The electron-vortex collision integral is diagonal 
when Fourier transformed, and we get the following equation for the n'th component of the 
distribution function 
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- evE^-^{5 nil + 6 n -i} = -n v {w - w n }f n - — — L^£. (52) 
The current is given by 

r fih p2^ R I Re[ ] 

(53) 




W 1 -W ) 

from which the conductivities can be read off. The resistivities are found by inverting the 
conductivity tensor. 

1. Longitudinal Resistivity 

The longitudinal resistivity, obtained by the above procedure, is 

m ( 1 1 \ 

Pxx = — [- + , (54) 

TIC y T v T~imp J 

where the transport scattering time for the electrons, due to scattering on the vortices, is 

1 = nvVF r dO{l - cos 6)\F(6)\ 2 . (55) 

Let us introduce a dimensionless quantity (, by writing the contribution to the longitudinal 
resistivity, from the electron-vortex scattering, as 

A. = (56) 

The dimensionless parameter (, which is explicitly given by the expression 

C = KjTd0(l-cos0)^, (57) 

is characterizing the efficiency of a single vortex, to scatter the electrons from the front to 
the back of the Fermi circle. It is straightforward to do the integral and obtain the following 
sum 



i=— 00 
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where t\ = tan 61 is found from Eqn. ^5|. 

In the Aharonov-Bohm limit we have an analytical expression for the cross section. It 
turns out that in this limit the integrand in Eqn. [55] is constant. We have in this limit 
C = 2 sin 2 vr/, or in other words 

1 sin 2 7r/ sin 2 71-/ 

— - = n v v F — = u c — , (59) 

t ab ttk 7i j 

where u c is the cyclotron frequency corresponding to the flux density of the external magnetic 
field. We note that ( vanishes in the Aharonov-Bohm limit, when the flux fraction / is 
integer, as it should. We can use this expression to make an estimate of the relative resistance 
change due to the vortices Ap xx /p xx = 2lffi v sin 2 Tif/kp. With / = 1/2, If = 5/im, n = 
10 n cm~ 2 and B = 100G, we get Ap xx /p xx = 0.6, which is a significantly higher value than 
is observed experimentally. Let us now include the effects of the finite radius of the vortices. 
Fig. ^| shows several ((k) curves for different values of the flux fraction /. It is seen that in 
general broader vortices give lower (, i.e. less resistance. The curve / = 1/2 corresponds to 
the physical Abrikosov vortices, if the difference in the cross sectional shape is ignored. In 
the very low field limit, that is less than 100G, the observed increase in resistivity is linear 
in the applied field, and the relative change at B = 100G is in the range 10~ 2 — 10~ 3 j3|, 
much less than the present estimate gives. 

The theory we have outlined here is only valid when the vortices are well separated, this 
amounts to assuming d 3> R v . The crossover to a different kind of behaviour observed in 
experiment, which is seen in both p xx and p xy , appear about 100G, and we believe that this 
is the flux density where the broad Abrikosov vortices begin to interfere. 



2. Transverse Resistivity 



The Hall resistivity, obtained from Eqn. IS3I, we can write as 



B , . 

Pxy = « — , (60) 

ne 
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where B is the externally applied homogeneous magnetic field. Here B/ne is the Hall 
resistivity of a 2DEG in a homogeneous magnetic field, and a is a dimensionless number, 
which describe the effects of the field being inhomogeneous. The dimensionless quantity a 
is given by the expression 

a = p- r d9 sin#|F(#)| 2 - 7^7 T dB sinfl^-. (61) 



Again, the integral can be expressed as a sum over terms involving only the parameters 
ti = tan 8i, which are given by Eqn. [45] 

Curves showing a as a function of k for different values of the flux fraction /, have been 
plotted in Fig. [8[ The classical limit is realized when R v /\p ^> 1 together with / ^> 1, and 
we have k = 27r-P 1L so none of the curves shown reach the classical regime. The number of 
terms, which must be included in the sum over /, grows rapidly with increasing k and /, 
thus making it difficult to reach the true classical regime by this technique. 

The a-curve for / = 1/2 we can compare with the experimental Hall factor measured 
by Geim et al. [[|. The overall qualitative behaviour is in good agreement, when the very 
idealized shape of the vortices, we have used in our calculation, is taken into account. To 
make a quantitative test we have fitted our curve to the experimental curve by tuning the 
radius of the vortex R v . The best fit is obtained for a vortex radius R v = 30nm, and this is 
nearly an order of magnitude smaller than the exponential length estimated by Geim to be 
lOOnm. This may indicate that most of the flux in the vortices are concentrated in a narrow 
core. 

It is seen in Fig. [8] that the a curves corresponding respectively to / = 1/4, and / = 
3/4 does not seem to converge to the level a = 1 as one would expect from the classical 
calculation. This is another manifestation of the Aharonov-Bohm phenomenon. Consider 
the ordinary double slit experiment with an infinitely thin solenoid hidden behind the middle 
obstacle. The interference pattern which can be observed behind the arrangement when it is 
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hit by an incident plane wave, is symmetric when the flux through the AB-solenoid is equal 
to zero, modulo the flux quantum. In this case the interference pattern will have a local 
maxima right in the middle. When the flux is equal to half a flux quantum, modulo the flux 
quantum, the interference pattern will again be symmetric, but this time with a node in 
the middle. In both the cases <fi = 0o an d 4> = 0o/2 there are no net scattering of electrons 
to either side. But for general fluxes f(f>o, with 2/ not equal to an integer, the interference 
pattern is not symmetric, and this is the reason that the / = 1/4 and / = 3/4 a-curves 
does not converge to the "classical" level at a = 1. When k ^> 1 and the scattering inside 
the vortex has become classically behaived, the vortex still plays the role of an obstacle that 
gives rise to an interference pattern, and the vector potential outside the vortex deflects the 
pattern and thereby gives rise to the asymmetry we observe in Fig. We emphasize that 
in the limit where the diameter of the vortex (the obstacle) goes to zero, there is no net 
scattering to either side. 

D. Multi Flux Quantum Vortex and Resonance Scattering 

When the total amount of magnetic flux inside the vortex is increased, the a and £ 
spectra acquire more structure. In Fig. |9] we have shown a and ( curves for a flux vortex 
carrying a total of 10 flux quanta. The structure seen in the plots is an effect of the resonant 
scattering which takes place when the energy of the incoming particle is close to one of the 
Landau quantization energies corresponding to the magnetic field strength inside the vortex. 
The magnetic field in the vortex vanishes outside a finite range - the radius of the vortex - 
and there are therefore no real Landau levels in the sense of stationary eigenstates, but only 
metastable states. In the dimensionless units we are working with, the Landau quantization 
energies E p = Tiuj c {p +1/2) corresponds to 

« P = 2^|/|(p+i), p = 0,l,2,... (63) 

These values are in excellent agreement with the resonances seen in Fig. [5], where the first 
eight resonances corresponding to p — 0, . . . , 7, are clearly distinguished. At the resonance 
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energies the typical time the particle spends in the scattering region, i.e. inside the vortex, 
is much longer than it is away from the resonance. The time the particle spends inside the 
vortex at a resonance, can be thought of as the lifetime of the corresponding metastable 
state. The inverse lifetime is proportional to the width of the resonance, that is strictly 
speaking the width of the peak in the partial wave cross section 07, corresponding to the I 
quantum number of the metastable Landau state. 

It is easy to interpret the small peaks in the £-curve in Fig. appearing at the resonance 
energies. Because when the electron spends longer time in the scattering region, it loses 
knowledge of where it came from, resulting in an enhanced probability of being scattered 
in the backwards direction. The a-curve is a measure of asymmetric scattering, and we 
can therefore interpret the dips seen in Fig. ^| along the same line of reasoning as for the 
C-peaks. The electron spends longer time in the scattering region, thereby losing knowledge 
of what is left and what is right. To explain why the a-curve is asymmetric in k around the 
resonances one could look at it this way: For increasing k the electron scattering becomes 
more and more classical, giving rise to the overall increasing background in a. But every 
time a new scattering channel is opened, the asymmetry of the scattering is suppressed, due 
to the lack of knowledge effect, thus resulting in a sawtooth like curve. 

III. HALL EFFECT IN A REGULAR ARRAY OF FLUX VORTICES 

A. Introduction 

Recently measurements were made by Geim et al. |2|,[|| of the Hall resistivity of low 
density 2DEG's in a random distribution of flux vortices, at very low magnetic field strengths. 
A profound suppression of the Hall resistivity was found, for 2DEG's with Fermi wavelengths 
of the same order of magnitude as the diameter of the flux vortices. This indicates that we 
are dealing with a phenomenon of quantum nature. These measurements were made by 
placing a thin lead film on top of a GaAs/GaAlAs heterostructure. When a perpendicular 



22 



magnetic field is applied, the magnetic field penetrates the superconducting lead film and 
also the heterostructure in the form of flux vortices each carrying half a flux quantum 0o/2 
of magnetic flux. Due to the strong flux vortex pinning in the films Geim have used, the 
vortices were positioned in a random configuration. 

In this section we will consider the hypothetical experiment where one instead of a "dirty" 
film, places a perfectly homogeneous type II superconducting film, on top of the 2DEG. The 
film do not have to be made of a material which is type II superconducting in bulk form. 
A film of a type I superconducting material will also display a mixed state if the thickness 
of the film is below the critical thickness d c . Experimentally perfect Abrikosov flux vortex 
lattices have been observed in thin films of lead with thickness d < d c = 0.1/im, [[17]]. 

If one succeeds to make such a sandwich construction, one has an ideal system for 
investigating how a 2 dimensional electron gas behaves in a periodic magnetic field. When 
the magnetic field exeeds H&, which can be extremely low, the superconducting film will 
enter the mixed phase, and form an Abrikosov lattice of flux vortices. The Abrikosov lattice 
in the superconducter will give rise to a periodic magnetic field at the 2DEG, and moreover 
as the strength of the applied magnetic field is varied the only difference at the 2DEG, is that 
the lattice constant of the periodic magnetic field varies. The Abrikosov lattice is most often 
a triangular lattice with hexagonal symmetry, although other lattices have been observed 
(e.g. square) in special cases where the atomic lattice structure impose a symmetry on the 
flux lattice, |[17|| . For simplicity the model calculations which we have done, were made for 



a system with a square lattice of flux vortices, but we do not expect this to influence the 
overall features of the results. From the point of view of the 2DEG it is important that the 
flux vortices carry half a flux quantum 4f = j- due to the 2e-charge of the Cooper pairs 
in the superconductor. The magnetic field from a single flux vortex fall of exponentially 
with the distance from the center of the vortex. This exponential decay is characterised by 
a length A s , which essentially is the London length of the superconducter, proportional to 
one over the square root of the density of Cooper pairs. The length A s can be varied by 
changing the temperature, or the material of the superconductor. 
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The other characteristic lengths of the system are the Fermi wavelength A^ = y ^f, where 
n is the density of the 2DEG, the lattice constant a of the periodic magnetic field, and the 
mean free path If = vft. The mean free path we assume to be very large compared to a 
and A s . The magnetic field is only varying appreciably when a is larger than A s . This means 
that the magnetic flux density of the applied field should be appreciably less than O / i^^l), 
which typically is of order lOOOGauss. In the limit where A s <C a, Xp, the vortices can be 
considered magnetic strings, and the electrons experiences a periodic array of Aharonov- 
Bohm scatterers. In this case the value of the flux through each vortex is crucial. If for 
instance the flux had been one flux quantum O — \, the electrons would not have been 
able to feel the vortices at all. But in the real world the vortices from the superconductor 
carry 4p = j 1 of flux, and therefore this limit is nontrivial. The electrons has for instance 
a band structure quite different from that of free electrons. In the mathematical limit of 
infinitely thin vortices each carrying half a flux quantum, there cannot be any Hall effect. 
This is most easily seen by subtracting one flux quantum from each vortex to obtain a flux 
equal to minus half a flux quantum through each vortex. As we have discussed earlier the 
introduction of the Dirac strings can not change any physics, and the procedure therefore 
shows that the system is equivalent to it's time reversed counterpart, thereby eliminating 
the possibility of a Hall effect. 

In this study we have ignored the electron spin throughout, in order to keep the model 
simple. From the point of view of the phenomena we are going to describe, the effect of the 
electron spin will be to add various small corrections. 

B. Electrons in a periodic magnetic field 

1. Magnetic translations 

It is a general result for a charged particle in a spatially periodic magnetic field B(x,y), 
that the eigenstates of the system can be labeled by Bloch vectors taken from a Brillouin 
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zone, if and only if the flux through the unit cell of the magnetic field is a rational number 
p/q times the flux quantum. The standard argument for this fact is made by introducing 
magnetic translation operators. To introduce magnetic translation operators in an inho- 
mogeneous magnetic field we first make the following observation. The periodicity of the 
field can be stated B(r + R) = B(r), for R belonging to a Bravais lattice. But this im- 
plies that the difference between the vector potentials A(r + R) and A(r) must be a gauge 
transformation 

V x {A(r + R) — A(r)} = B(r + R) - B(r) = 0. (64) 

We introduce the gauge potential xr an d write 

A(r + R) = A(r) + V XR (r). (65) 

The function xr is only defined modulo an arbitrary additive constant which have no physical 
effect. The Hamiltonian of the electrons is 

H = ^-(p + eA(r)f. (66) 

The ordinary translation operators Tr = exp[^R-p] do not commute with the Hamiltonian, 
because they shift the argument of the vector potential from r to r + R, but as we just have 
seen this can be undone with a gauge transformation. We therefore introduce the magnetic 
translation operators as the combined symmetry operation of an ordinary translation and a 
gauge transformation 

6 i 

M R = exp[-i-Xfl(r)] exp[-R ■ p\. (67) 

The operator Mr is unitary, as it is the product of two unitary operators, and therefore has 
eigenvalues of the form e lX . Let us denote the primitive vectors of the Bravais lattice a and 
b. We can find common eigenstates of M a , M b and H, if and only if they all commute with 
each other. The magnetic translations each commute with the Hamiltonian by construction, 
and furthermore we have 
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M a M b = exp[2%i-^}M b M a . (68) 
<Po 

If the flux through the unit cell is a rational number p/q (p and q relatively prime) times 
the flux quantum <p Q , M qa and M b commute. In this case the cell spanned by qa and b 
is called the magnetic unit cell. Let us define c = qa. The possible eigenvalues of M c are 
phases e 2mkl , where we can restrict \ki\ < |, and equivalently for M b . We can therefore label 
the common eigenstates \k,n), where k = k\C* + k 2 b*, and c*, b* are the primitive vectors of 
the reciprocal lattice. The vector k is restricted to the magnetic Brillouin zone. An arbitrary 
magnetic translation of an eigenstate with a Bravais lattice vector R = nc + mb can now be 
written M R \k,n) = (M c ) n (M b ) m \k, n) = exp[ik ■ R]\k,n), showing that the eigenstate is a 
Bloch state. In this case we can speak of energy bands forming a band structure in the usual 
sense. When the flux through the elementary unit cell (a, b) is an irrational number of flux 
quanta, the situation is different. The irrational number can be reached as the limit where 
p and q get very large, and consequently the Brillouin zone get very small and collapses in 
the limit. 

2. The Dirac vortex viewpoint 

In this section we will show, how it is possible to argue in a slightly different way from 
the previous section, and hereby in a simpler way obtain the vector potential of a periodic 
magnetic field. Let us again assume a rectangular unit cell (a, b), B(x + a,y) = B(x,y + 
b) = B(x,y) etc., to keep the notation simple. The magnetic field enter the Hamiltonian 
only through the vector potential. The question is therefore if one can choose a gauge 
such that the vector potential will be translationnally invariant relative to a unit cell (c, d) 
A(x + c, y) = A(x,y + d) = A(x, y) etc. It is clear that if such a periodic A-field exists, then 
the total flux <3> crf through the unit cell (c, d) will be zero, as it is given by the line integral 
of A around the boundary of the unit cell (c, d) 

= / A dl, (69) 
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which is zero by the periodicity. We remark that due to the relation B = V x A, the cell 
(c, d) will be bigger than or equal to (a, b). If the flux through the unit cell of the magnetic 
field is not zero, but equal to a rational number times the flux quantum, $ a f, = p/q (f> [p 
and q relatively prime), a trick can be applied to make the flux $ c d become zero. It is a 
basic fact, apparently first observed by Dirac [|18|| , that a particle with charge e cannot feel 
an infinitely thin solenoid carrying a flux equal to an integer multiple of the flux quantum 
0o — h/e. Such a stringlike object carrying one flux quantum is sometimes called a Dirac 
vortex. On a lattice the Dirac vortex goes through the center of a plaquette, and the electron 
can therefore never enter the core of it. The lattice Dirac vortex can be made to disappear 
by a gauge transformation. 

To find the periodic vector potential take an enlarged unit cell (c,d) = (qa,b), so that 
&cd — P<Po and put by hand p counter Dirac vortices through the cell, to obtain zero net 
flux. Then a divergence free vector potential can be build for instance by Fourier transform 



B(Q) = — / d 2 rexp[-iQ ■ r)B(r), (70) 

Cd J(c,d) 

B(Q), (71) 



Mr) = £ 



(c,d) 

\ 

exp[iQ ■ r] 



Q' 1 



\ -iQx j 

where the sum is over Q in the reciprocal lattice. Here we have used continuum notation, 
but it is straightforward to write down the lattice equivalents of the expressions. 



C. Lattice calculation of Hall conductivity 

We have calculated the Hall conductance of the 2DEG in the vortex field by a numerical 
lattice method. This we do because the calculations then reduces to linear algebra operations 
on finite size matrices, which can be implemented in a C++ program on a computer. The 
idea is to consider an electron moving on a discrete lattice, rather than in continuum space. 
We know, allthough we are not going to prove it here, that in the limit where the discrete 
lattice becomes finegrained compared to all other characteristic length of the system, the 
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continuum theory is recovered. Here we assume that the original Bravais lattice has square 
lattice symmetry, with a lattice parameter which we call a. The discrete micro lattice is 
then introduced as a finegrained square lattice inside the unit cell of the Bravais lattice. The 
lattice parameter of the micro lattice we then take as a/d, where d is some large number, in 
order to keep the to lattices commensurable. The condition that we have to impose on the 
micro lattice, in order that it is a good aproximation to the continuum, can then be stated 

a/d <C a, Xf, A s , • • • . (72) 

In the numerical calculations we have made, we have taken d = 10. 
The tight-binding calculations are made with the Hamiltonian 

H = — ^2 ti+T,j+T'C i+T Cj +T >. (73) 

ijrr f 

Here i, j are Bravais lattice vectors, and r, r' are vectors indicating the sites in the basis. 
The matrix elements t, +r j +T > are taken non-zero only between nearest neighbour sites. The 
matrix element between two nearest neighbour sites r and t' are complex variables t TtT + e = 
j-giA^r) w j^} 1 a phase given by the vector potential A^{r) residing on the link joining the sites. 
The translation invariance of the Hamiltonian can then be stated ti +T+ ij +T i + i = ti +T j +T i, 
for all vectors I belonging to the Bravais lattice. Let us introduce the system on which our 



calculations were made as an example. Fig. 10 shows the unit cell with its internal structure 



i.e. the basis. There are N = d ■ d sites in the basis. The length of the links we write 
as a/d, where a is the side of the unit cell, with area Q = a 2 . The vectors r = (ti,T2)§, 
7~i,T2 = 0,1,... d— 1 are offsets into the basis, while the vectors i = {i\,i2)a, € Z 
indicate the cells in the Bravais lattice. The operator cj +T , for a given r, is defined on the 
Bravais lattice, and accordingly it can be resolved as a Fourier integral over the Brillouin 



4. = 7tW ( 74 ) 



zone as 

Q 

(2^) 2 Jbz' 

(It should be noted that the factor e~ tq " T is arbitrary and included here for later convenience). 
Inserting this and using the translation invariance, the Hamiltonian can be rewritten as 
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H =w?L d2kH *' < 75) 

where we have introduced 

= " t T , +T/ e lk -^'-^cl T c ky . (76) 

It is seen that Hk only mixes the N states \kr), i.e. it is an A^ x A" matrix. The A" eigenvalues 
of Hk are the energies of the A" tight binding Bloch states with wavevector k. Let us denote 
the eigenstates of Hk by 

H h ul = Etui (77) 

where a — 1,2, ... N and E% < E^ +1 - From the A" dimensional vector we can construct 
the eigenstate of the Hamiltonian H 

(i + r|^) = e ife -^<(r). (78) 

It is straightforward to verify that this is the correct Bloch eigenstate of H. The bandstruc- 
ture can be calculated directly by diagonalising the N x N matrices, Hk, for representative 
choices of k in the Brillouin zone. Before one can compare the spectrum optained from 
this calculation with that of a continuum system, a scaling of the energies is required. To 
scale the energy to the spectrum of a particle with an effective mass m, we have to take 
t = h 2 d 2 /ma 2 , and e (fe) = E% + At. 

The Hall conductivity can be calculated by the same method as in the homogeneous 
magnetic field |7|,[T9|. We have used a single particle Kubo formula to calculate the Hall 
conductance 



a 



ifl \ - {Jx)ap{Jy)l3a ~ (Jy) a p{Jx) 0a / ?q x 



' A ° (E« - BP) 



where J x ,J y are the currents in the x,y directions, and the sum is over single particle states 
\a, k) with energies below and above the Fermi level Ep. The area of the system is denoted 
Aq. All quantities are diagonal in k, and therefore this index is suppressed. The summation 
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is composed of a discrete sum over bands, and an integral over the Brillouin zone for each 



band. The Brillouin zone shown in Fig. |11| is doubly connected because the states on the 
edges is to be identified according to the translation invariance. This gives the Brillouin zone 
the topology of a torus T 2 , with two basic non contractible loops. The current operator can 
be written 

where Jk = § ^f- By use of some simple manipulations and completeness, it is straightfor- 
ward to refrase Eqn. [79| 



where 



162 - (lda - (81) 



da \ I da 



da 



dk x / / ' 

■^-^j is shorthand for g|- \(3,k). This formula was first derived by Thouless, Kohmoto, 
Nightingale and den Nijs [[?]], for a noninteracting 2 dimensional electron gas in a periodic 
scalar potential, and a commensurate perpendicular magnetic field. It requires some com- 
ments to be meaningful. In order to calculate -ir-) it is necessary to consider the difference 



(\a, k + 5k^) — | a, k))/Sk fl . But this difference is not well defined as it stands, as the phase 
of the states is arbitrary. Rather than representing the state by a single vector in C N , it 
should be represented by a class of vectors which differ one from another only by a phase. 
These equivalence classes are sometimes called rays. To compare states locally, we need to 
project this U(l) degree of freedom out. This is done by demanding the wave function to 
be real, when evaluated in a fixed point, i.e. u^{t{) = (xj|o;,fc) e R. If the wave function 
happens to be zero in r iy some other point Tj must be used. When a band has a non- zero 
Hall conductivity, it is not possible to find a single r which work for all the states in the 
Brillouin zone. The change from Tj to Tj which shifts the phase of the states, is analogous to 
a gauge transformation on the set of states. The special combination of terms which appear 
in Eqn. [79| is gauge invariant with respect to these special "gauge transformations" . If we 
let \x a ) denote a state which is obtained from \a,k) by fixing the phase according to the 
above scheme, the following formula for the contribution to the Hall conductivity from a 
single band a, is well defined 
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dk„ 



dx c 



dk x 



dk x 



d X a 



dk„ 



(82) 



It has been shown in detail by Kohmoto [pl|], for the homogeneous magnetic field case, 

2 

that this expression is equal to minus %- times the first Chern number of a principal fibre 
bundle over the torus. As the first Chern number is always an integer, this has the physical 
consequence that whenever the Fermi energy lies in an energy gap, the Hall conductance 
is quantized. We will use this result to interpret certain peaks in the cx^-spectra we have 
calculated. When the Fermi level is not in an energy gap of the system, we will have to use 
Eqn. |8l] to calculate a xy . As we shall see, in this case there is no topological quantization of 
the Hall conductivity. 



D. Energy band crossing 

In this section we study the effect on the Hall conductivity of an energy band crossing. 
This has previously been discussed in different contexts by several authors [|20| |22 1 . 



When the shape of the magnetic field is varied, controlled by some outer parameter £, 
it will happen for certain parameter values £o> that two bands cross, see Fig. IT2]. This is 



the consequense of the Wigner-von Neumann theorem, which states that three parameters 
are required in the Hamiltonian in order to produce a degeneracy not related to symmetry. 
Here the parameters are k x , k y and the outer parameter £, which in our calculation is the 
exponential length of the flux vortices from the superconducter. When the energy difference 
E + — E~ between the two bands considered is much smaller than the energy distance to the 
other bands, the Hamiltonian can be restricted to the subspace spanned by the two states 
and | — , A; ). The point in the Brillouin zone where the degeneracy occur we denote 
k°. The Hamiltomian H(k) is diagonal for k = k°, and we denote the diagonal elements 
respectively E + e and E — e. For small deviations of k from k° the lowest order corrections 
to the Hamiltonian is offdiagonal elements A(fc) linear in k — k°. Without essential loss 
of generality we can assume that e is independent of k. Then e plays the role of the outer 
parameter controlling the band crossing. The Hamiltonian is then approximated by 
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H(k) = 

The offdiagonal element is expanded as 



e A* 



A -e 



+ E . 



(83) 



A(fe) = a{k x - k° x ) + P(k y - k° y ) 
with a = k Q \H k \+, fc°>, and /3 = k°\H k \+, k°). 



(84) 



We want to find the consequenses of the energy band degeneracy, on the topological Hall 
quantum numbers of the bands. Let us define 



B ± (k) 



'd± 
dk v 



d± 



d± 



d± 



dk„ 
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Then the interesting quantities are the integrals of B + (k) and B-(k), around a small neigh- 
bourhood of the degeneracy point fe°. It turns out that it is the two numbers a and (3 that 
control what happens. 

In general a and f3 will be nonzero complex numbers — nonzero because we have assumed 
the degeneracy to be of first order. Let us first consider the degenerate case where a and (3 
are linearly dependent, i.e. a/ (3 is real, or otherwise stated lm(a*/3) = 0. Then by a linear 
transformation we can write the Hamiltonian 



(86) 



h(K) = (ki + K 2 )cr 1 + 70- 3 = 

^ K X + K 2 -7 J 

where now 7 is the dimensionless parameter of the crossing, and k±,k 2 are the rescaled 
dimensionless momentum variables. The <r M 's refers to the Pauli matrices. This Hamiltonian 
is real, and we can therefore also choose the eigenstates to be real, and this will clearly lead to 
a vanishing B(n). In this case we therefore conclude that there is no exchange of topological 
charge. Here the word topological charge is used to denote Hall quanta. 

Let us now treat the general case where a and f3 are linearly independent, i.e. lm(a*fl) 7^ 
0. In this case, we can by a linear transformation write A(k)/E — k± + in 2 = ne l6 ', which 
defines the scaled momentum variables k, 9. (Here Eq is some constant with dimension of 
energy.) This reduces the Hamiltonian to the form 
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h(n, 0) = KiiT 1 + k 2 <7 2 + l^ 3 



' 7 ne-"^ 



K,e 



w 
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-7 J 



Let us define A = a/7 2 + k 1 . Then the eigenvalues of h(n, 9) are ±A and the two correspond- 
ing eigenstates are 



/ 



|±,fc) 



y/2 



±J1± 7 /A 
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V 
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In order to calculate the integral of the 5-function we need to express it in terms of the 
k, variables. The Jacobian of the transformation is given by the expression 

k dndO 



d k% d ky 



(89) 



and 
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d± 
dn 



d± 
dn 



^)} = ±KA-«^)^, (90) 



where the indices r, i refer to the real and imaginary parts respectively. We can now calculate 
the contribution to the Hall conductivity from each of the bands, from the area around k° 
given by \k\ < k c , where k c is some local cutoff parameter which limit the integration to the 
area where the approximation leading to the Hamiltonian Eqn. |83| is valid 
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where the last equality sign is valid when n c /\l\ ^> 1, i-e. close to the crossing where 
7 = 0. Here the factor sign [7] signals that the two bands exchange exactly one topological 

2 

conductivity quantum at the crossing. This is not surprising, because we know that the 
total contribution from the states in a single band is always an integer times the quantum 
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Moreover it is readily seen that in the hypothetical situation of a n'th order degeneracy, 
i.e. one for which A = K, n e me , n quanta are exchanged. The total topological charge in 
the band structure is conserved. The topological charge can flow around and rearrange 
itself inside a band, but only be exchanged between bands in lumps equal to an integer 

2 

multiplum of the conductivity quantum j-. When the radius of the flux vortices is gradually 
shrinked to zero, the Hall effect has to disappear. There are two mechanisms with which 
the Hall effect can be eliminated. The first is by moving the topological charge up through 
the band structure by exchanging quanta, resulting in a net upward current of topological 
charge, eventually moving the charge up above the Fermi surface, where it has no effect. 
The second mechanism is by rearranging the topological charge inside the bands, so that 
each band has a large negative charge in the bottom, and a large positive charge in the 
top, but arranged in such a clever way that charge neutrality is more or less retained for 
all energies. This second mechanism will also give a net displacement of topological charge 
up above the Fermi energi, because in general the Fermi surface cuts a great many bands, 
and for all these bands the large negative charge, which they have in their bottom part, will 
be uncompensated. To use the language of electricity theory, we can say that every band 
gets extremely polarized, resulting in a net upward displacement current, in analogy with 
the situation in a strongly polarized dielectric. Our numerical calculations indicate, that it 
is the second mechanism which is responsible for the elimination of the Hall effect, as the 
radius of the vortices shrinks to zero. 

When two bands are nearly degenerate for some k°, each of the bands have concentrated 
half a quantum in a small area in k-space around k°, and in general the topological charge 
piles up across local and global gaps in the energy spectrum. This is the reason for the 
oscillatory and spiky behaviour of the Hall conductance as a function of electron density, that 
is seen on the calculated spectra below. It is also the reason why the numerical integration 
involved in the actual evaluation of the Hall conductivity, is more tricky than one could 



wish. In particular it indicates that it is not true as sometimes conjectured, e.g. [23], that 
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the Hall conductivity is smoothly distributed in k-space, and that it therefore should vary 
smoothly between the quantized values at the energy gaps, as the Fermi energy is swept 
through a band. 



E. Numerical results 



1. Transverse Conductivity 



We have calculated the transverse conductivity (7 X y cLS cL function of the integrated density 
of states, for electrons in a square lattice of flux vortices, for a series of varying cross 
sectional shapes of the flux vortices. Each of the field configurations consists of a square 
lattice of flux vortices with a given exponential length A s . The parameter which vary from 
calculation to calculation, is the dimensionless ratio £ = X s /a, where a is the length of 



the edges of the quadratic unit cell. The unit cell is shown on Fig. |10| and contain, as 
we have already discussed, two vortices and a counter Dirac vortex. In order to do the 
tight binding calculation a micro lattice is introduced in the unit cell. In all the numerical 
calculations we present, the micro lattice is 10 x 10. This gives 100 energy bands distributed 
symmetrically about the center on an energy scale. Out of these only the lower part, 
say band 1 to 20, approximate the real energy bands well, while the rest is significantly 
affected by the finite size of the micro lattice. A careful examination of the vector potential 
reveal that the symmetry of the Hamiltonian is very high for the particular choice of unit 
cell shown in Fig. [H| The field from a single vortex we have taken as Boe^^ Tx ^ + ^ Ty ^^ Xs 
instead of the more realistic Boe~^^ Xa . With this choice the vector potential can be written 
down analytically in closed form. This makes the calculations simpler, and does not break 
any symmetry that is not already broken by the introduction of the micro lattice. (We 
have made calculations of the band structure, with both kinds of flux vortices, and the 
differences are indeed very small). The energy spectrum is invariant under the changes 
(k x , k y ) i — > (±k x , ±k y ), {±k y , ±k x ). This fact is exploited to present the band structures in 
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an economic way. The labels T, X and M correspond to the indicated points in the Brillouin 



zone, Fig. 11 



In Fig. ^ we have plotted a selection of typical bandstructures which illustrates the 
crossover from the completely flat Landau bands in the homogeneous magnetic field £ = oo, 
to the bandstructure of electrons in a square lattice of Aharonov-Bohm scatterers (/ = 1/2) 
at £ = 0. The bandstructures have been found by direct numerical diagonalization of the 
Hamiltonian. (See also Fig. |16|). 

In Fig. [14] some typical results of the numerical calculations of the Hall conductivity 
are shown. In general we have no reason to expect, that the Hall conductivity should be 
isotropic as a function of the angle between the current and the flux vortex lattice. The 
results we present is for a current running along the diagonal of the square lattice, i.e. along 
the x-axes in Fig. |10[ It is seen, that whenever there is a gap in the spectrum, the Hall 
conductivity gives the quantized value in agreement with the discussion in the last section. 
At Fermi energies not lying in a gap o xy allways tend to be lower than the value it has in 
the homogeneous field. And in the limit £ i — > 0, a xy vanishes alltogether. In this limit 
the electron sees a periodic array of Aharonov-Bohm scatterers, each carrying half a flux 
quantum, and there is no preferential scattering to either side. We observe that when the 
flat Landau bands starts to get dispersion, the contribution to the Hall effect is no longer 
distributed equally in the Brillouin zone. Instead it piles up across local and global gaps in 
the spectrum resulting in the spiky a xy spectra. Finally, we note that the amplitude of the 



2 



fluctuations of a xy is of order j-, consistent with this picture. 



In the calculations presented in Fig. |14], the filling fraction is limited to values below 
30. This limitation comes from the fact, that we are only able to handle matrices of limited 
size in the numerical calculations. Filling fractions below 30 correspond to very low density 
electron gases, and our calculations belong therefore to the "quantum" regime, i.e. to the 
regime where A^ 3> X s . According to the discussion in Sec. [TT], we expect a cross-over 
to a semiclassical regime, for electron gases with higher density, where A^ <C A s , with a 
qualitatively different behaviour. 
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It is an important question wether it is possible to observe these features of the Hall 
conductivity in experiments. The conditions, which are nescessary, are that the mean free 
path I is long compared to all other lengths, and that a, X s < Xp. The last condition 
indicates that we are concerned with quantum magneto transport, in the sense that the 
vector potential is included in the proper quantum treatment of the electrons. This is 
in contrast to the case -C a, X s where the electron transport can be treated as the 
semiclassical motion of localized wavepackets in a slowly varying magnetic field. Let us 
assume, in order to make some estimates, that the superconductor has a London length 
somewhat less than 1000A, resulting in an exponential length of the vortices X s of about 
1000A at the 2DEG, after the broadening due to the distance between the superconducter 
and the 2DEG has been taken into account. In order to have a variation in the magnetic 
field we should have a > X s , and to be in the quantum regime Xp > a, X s . This gives the 
estimate for the electron density n < 6 • 10 10 cm -2 , which is not unrealistic. The effect of 
the impurities is (to first order), to give the electrons a finite lifetime. This gives a finite 
longitudinal conductivity a xx and broadens the density of states. It also introduces localized 
states at the band edges (Lifshitz tails). If the field is nearly homogeneous, we have the 
standard quantum Hall picture with mobility edges above and below every Landau band, 
resulting in the formation of plateaus in a xy , which only can be observed at much higher 
magnetic fields, of order 10 5 G, where the filling fraction is of order one. On the other hand 
when the amount of impurities is low, that is kpl 3> 1, all these effects will be small, and 



we expect that the features of a xy , shown in Fig. 14 will have observable consequenses in 



Px y °i*+oiy' 



In Fig. 16 we have plotted the bandstructure of electrons in a square lattice of vortices 



carrying one flux quantum each. The calculation has been made with the same basis as the 



bandstructures shown in Fig. [13], the only difference is that all fluxes have been multiplied 
by a factor of 2, as can be seen from the spacing between the Landau bands. In the limit 
where £ — > 0, and the vortex lattice becomes a regular array of Aharonov-Bohm scatterers, 
we recover the well-known bandstructure of free electrons. This is in agreement with our 
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discussion of the AB- vortex in Sec. ||. 



2. Exchange of topological quanta 

An example of exchange of topological quanta between neighbouring bands is shown in 
Fig. [T7|. The figure is an enlargement of the band structure around the X point, showing an 
accidential degeneracy between the 3'ed and 4'th band, which occur about £ = 0.035. Also 

2 

indicated is the Hall conductance of each band in units of found by numerical integration. 
At the degeneracy it is not possible to define the Hall conductance for the individual bands. 
On the Brillouin zone torus there are two X points X%, X2, and the bands have a simple (1. 
order) degeneracy in each. The numerical integration shows that two topological quanta are 
transfered from the lower to the upper band, and this is in full agreement with the discussion 
of Sec. |H D| . Exchange of topological quanta between bands is a common phenomena as £ 



is varyed, and this particular example has only been chosen as an illustration of the general 
phenomena. 

3. Transverse conductivity in a disordered vortex phase 

When the vortices come from a thin film of superconducting material, which have many 
impurities and crystal lattice defects acting as pinning centers for the vortices, the distri- 
bution of vortices will be disordered rather than forming a regular Abrikosov lattice. The 
effect of the disorder will be to wash out the destinctive features of the band structure, 
i.e. to average out the characteristic fluctuations in a xy , leaving a smooth curve in Fig. |l4] 
with a characteristic dimensionless propertionality constant s(£), in the form a xy = ^-s(£)i/, 
where v is the number of electrons in the magnetic unit cell v = na 2 = (the filling 
fraction). With this conjecture we can estimate the normalized Hall conductivity s(£) by 
making a linear fit to the calculated c xy [y) distribution. In the experimental situation A s is 
constant, and this makes s(£) a function of the applied magnetic field through the relation- 



ship £ = X s /a = \ S J B/<j)Q. In Fig. 18 we have plotted s(B) for a vortex exponential length 

38 



X s = 80nm. The s(I?)-curve shows essentially that the Hall effect of a dilute distribution of 
vortices is strongly suppressed compared to the Hall effect of a homogeneous magnetic field 
with the same average strength. This is in good qualitative agreement with what is seen in 
the experiments of Geim et al. 0. When doing experiments, one is not directly measuring 
the conductivities, but rather the resistivities p xx , p xy . The experiments of Geim et al. cover 
the parameter range from A^ <C A s at high 2DEG densities, down to the value Xf/X s = 0.7 
for the 2DEG with the lowest density experimentally obtainable, where the new phenom- 
ena begin to occur. Our numerical calculations belongs to the other side of this cross-over 
where Xf ^> X s . The physical picture of this cross-over can be stated as follows. On the 
high density side the magnetic field varies slowly over the size of an electron wavepacket for 
electrons at the Fermi energy, with the result that the wavepacket more or less behaves as 
a classical particle. On the other side of the crossover Xp ^> X s the magnetic field varies 
rapidly over the lengthscale of a wavepacket, for an electron at the Fermi energy, and this 
introduces new phenomena of an essential quantum character. 



F. Superlattice potential 

The general picture we have outlined so far of energy bands having dispersion, with the 
dispersion giving rise to a non trivial behaviour of the Hall conductivity, is not limited to 
the inhomogeneous magnetic field. The dispersion could have another origin for instance a 
superlattice potential. To illustrate this a series of calculations have been made on a 2DEG 
in a homogeneous magnetic field, and a scalar potential which we have taken as a square 
lattice cosine potential. This system have commensurability problems because of the two 
"interfering" length scales, given respectively by the magnetic length Ib = and the 

period of the superlattice potential a. To make things as simple as possible we have fixed 
the period of the cosine potential a, and the magnetic field strength B, and only varied 
the amplitude of the cosine potential. Furthermore the flux density of the magnetic field 
is tuned so that the flux through one unit cell of the cosine potential is exactly one flux 
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quantum. The cosine potential is 

U(x,y) = Vo(cos27rx/a + cos2-7n//a), (92) 

and the magnetic field B = (po/a 2 . The dimensionless parameter controlling the shape of 
the energy band structure is in this case 

» = (93 > 

Examples of the cx^-spectra are shown in Fig. [19|. It is observed that although the spectra 
look different from the vortex lattice spectra, they have the same spiky nature. The spikes 
have the same interpretation as in the vortex lattice system. Local spikes are due to local 
gaps in the spectra. That is when to bands are close to each other for some k vector in the 
Brillouin zone, the result is a pile up of topological charge across the gap, and this gives a 
spike in the o^-spectra when the Fermi energy is swept across the gap. Global spikes, that 
is spikes which go all the way up to the diagonal line indicating the Landau limit, are due 
to global gaps in the energy spectrum, combined with the topological quantization. 



IV. CONCLUSION 

In Sec. H °f this paper we have considered the longitudinal and the transverese resistivi- 
ties of a 2DEG in a disordered distribution of flux vortices, within the theoretical framework 
where each scattering event is treated independently, and the electrons are non-interacting. 
The general features observed in experiments are in agreement with the results we have out- 
lined, but we do not have quantitative agreement. If we use the radius R v of the vortices as a 
fitting parameter, then the longitudinal resistance fits to a radius of the order R v ~ lOOOnm, 
while the transverse resistivity fits best at R v ~ 30nm. The radius of the real vortices, is 
estimated by Geim to be R v ~ lOOnm. This point requires analysis on a more elaborate 
level, in order to be resolved. 

In Sec. [Til] we have considered a new kind of experiment where a 2DEG is placed in a 
periodic magnetic field varying on a length scale A s , comparable to (or less than) the Fermi 
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wavelength of the electrons. In this limit, where it is nescessary to include explicitly the 
vector potential in a quantum treatment of the electron motion, we expect the 2DEG to 
exhibit new phenomena. We have presented numerical results for a non-interacting 2DEG 
without impurities showing characteristic spikes of the Hall conductivity versus filling frac- 
tion, which can be understood in terms of local and global energy gaps in the spectrum. 
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FIGURES 

FIG. 1. The scattering geometry for classical scattering on an idealized cylindrical vortex with 
constant magnetic field inside and no field outside. 

FIG. 2. Differential cross section and classical trajectories for four different values of the pa- 
rameter 7 = l c /R v = 0.025, 0.40, 1.00, 2.50 

FIG. 3. The geometry of the scattering situation. 

FIG. 4. The dimensionless Aharonov-Bohm cross section, defined as T AB (9) = J2i •7 7 [<5/ LB ]e de '. 
Here plotted for / = 1/2. The different curves are cross sections corresponding to different values 
of k, and they have been translated relative to each other in order not to overlap to much. The 
horisontal lines indicate the zero level for the different curves. The lowest curve corresponds to 
k = 0.25, and the other curves to respectively 1.00, 2.50, 5.00, with k increacing upwards. The 
scale of the plot can be read off the distance between the horisontal lines, which is equal to 1. 

FIG. 5. Differential cross sections plotted relative to the Aharonov-Bohm cross section, for an 
electron scattering on a magnetic flux vortex with finite radius. The curves have been obtained 
by subtracting the dimensionless AB cross section from the calculated finite radius cross sections. 
This have been done in order to single out the effect of the finite radius. All plots are in the same 
scale, and this includes the preceding figure, see caption of Fig. 4. Furthermore all plots are with 
the same values of n = 0.25, 1.00, 2.50, 5.00 increasing upwards. 

FIG. 6. Resistance efficiency of single vortex. Curves show ((k) for different values of the flux 
fractions /. 

FIG. 7. Resistance efficiency C(t)> an d Hah 1 efficiency factor a(j) calculated from the classical 
cross section for scattering on a magnetic flux tube. The parameter 7 is given by the cyclotron 
radius divided by the radius of the flux tube 7 = l c /R v . 
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FIG. 8. The efficiency of a dilute distribution of vortices in producing Hall effect, compared to 
a homogeneous magnetic field with the same average flux density. Curves show a as a function of 
re for different values of the flux /. 

FIG. 9. These plots of a and ( for a vortex with / = 10, shows a striking structure of resonances 
at the values of re corresponding to the Landau quantization energies. 

FIG. 10. (A) The unit cell with basis. The large circles indicates the position of the Abrikosov 
vortices, and the small circle indicate the position of the Dirac vortex with the counter flux. The 
micro lattice shown here is 6 x 6, whereas all the numerical results we have presented are obtained 
with a micro lattice of 10 x 10 sites. (B) Four concatenated unit cells, showing the squarre lattice 
of Abrikosov vortices. 

FIG. 11. The Brillouin zone with the conventionel symmetry labels. 

FIG. 12. Schematic energy band crossing, controlled by an outer parameter 7 = £ — £q- 

FIG. 13. Bandstructures for 2D electrons in a square lattice of Abrikosov vortices with / = 1/2. 
The different plots show bandstructures corresponding to various values of the parameter £ equal 
to the ratio between the exponential length of the magnetic field from a single vortex, and the 
lattice parameter. The flux through a single vortex is half a flux quantum. 

FIG. 14. Calculated Hall conductivity versus filling fraction, for various values of the ratio 
£ = \ s /a. These calculations are made on the same system as the bandstructures of Fig. 3.4, that 
is a square lattice of Abrikosov vortices with / = 1/2. Each of the spectra are made as follows. For 
2000 equidistant Fermi energies cf, the total Hall conductivity <jh(£f), and the integrated density 
of states v(€f) are calculated numerically. This is done by 20 pages of C++ code, running on a 
workstation for 24 hours. The x-axes indicates the integrated density of states in units of filled 
bands. The y-axes indicates the total Hall conductivity in units of the conductivity quantum e 2 /h. 
The diagonal line in the plots indicate the Hall conductivity in a homogeneous magnetic field. 
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FIG. 15. The density of Hall effect, or "topological charge", plotted as function of the filling 
fraction. It is seen that for £ — ► the distribution gets strongly polarized, with a negative con- 
tribution to the Hall effect at the bottom part of the bands, and a positive contribution at the 
uppermost part of the bands. 

FIG. 16. Bandstructures for 2D electrons in a square lattice of vortices carrying one flux quan- 
tum each, / = 1. The different plots shows bandstructures corresponding to various values of the 
parameter £ equal to the ratio between the exponential length of the magnetic field from a single 
vortex, and the lattice parameter. The bandstructures have been calculated using the basis shown 
in Fig. 10, with the only difference that here the flux through each of the vortices are fo, and the 
counter Dirac flux is — 2<^>o- 

FIG. 17. Exchange of topological quanta. The figure shows an enlargement of the / = 1/2 
band structure around the X point, where a degeneracy between the 3'rd and 4'th band occur. 
(See Fig. 13). The parameter 7 appearing in the figure is defined as 7 = £ — £o 5 with £0 = 0.035. 

2 

The numbers give the Hall conductance of the bands in units of found by numerical integration. 

FIG. 18. The normalized Hall conductivity s(B) for vortices of exponential length X s =80nm. 

FIG. 19. Calculated Hall conductivity versus filling fraction, for a 2DEG in a homogeneous 
magnetic field, and a square lattice cosine potential, in the special case where the magnetic flux 
density is exactly equal to one flux quantum per unit cell area. The dimensionless parameter 
v indicated in the plots, is equal to the amplitude of the cosine potential divided by the Landau 
energy hu c . The x-axes indicates the integrated density of states in units of filled bands. The y-axes 
indicates the total Hall conductivity in units of the conductivity quantum e 2 /h. The diagonal lines 
in the plots indicate the Hall conductivity in a homogeneous magnetic field, without any potential. 
For further details see the caption of Fig. 14 
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